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Abstract 



We consider several renormalizable, scale free models in three space-time 
dimensions which involve scalar and spinor fields. The Yukawa couplings are 
bilinear in both the spinor and scalar fields and the potential is of sixth or- 
der in the scalar field. In a model with a single scalar field and a complex 
Fermion field in three Euclidean dimensions, the couplings in the theory are 
both asymptotically free. This property is not retained in 2 + 1 dimensional 
Minkowski space, as we illustrate by considering a renormalizable scale-free 
super symmetric model. This is on account of the different properties of the 
Dirac matrices in Euclidean and Minkowski space. We also examine a model 
in 2 + 1 dimensional Minkowski space in which two species of Fermions, asso- 
ciated with the two unitarily inequivalent representations of the 2x2 Dirac 
matrices, couple in two different ways to two distinct scalar fields. There are 
two types of Yukawa couplings in this model, and either one or the other of 
them can be asymptotically free (but not both simultaneously). 

I. INTRODUCTION 

Self-interacting scalar models are the simplest of all relativistic quantum theories to 
analyze. There are but three instances in which the self-interaction involves a dimensionless 
coupling A, (and is consequently perturbatively renormalizable): in four dimensions, 

in six dimensions and A0^ in three dimensions. The renormalization group functions for 
these models have been determined beyond one-loop order using dimensional regularization 

see refs. respectively. In this paper we consider not just sixth order scalar couplings 
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in three dimensions, but append to them Yukawa couphngs which are bihnear in both the 
scalar and spinor fields. This interaction is renormalizable in three dimensions. 

We compute the renormalization group functions associated with these couplings to first 
order in the scalar self-coupling and second order in the Yukawa coupling in a variety of 
models. First, a simple model in which there is a single scalar field and a two-component 
Dirac Fermion field ip is considered in Euclidean space and we find it to be asymptotically 
free. Next a supersymmetric model in Minkowski space is examined and, on account of 
the properties of Dirac matrices in this space, the model no longer possesses asymptotic 
freedom. Finally, a model containing two species of two-component Dirac Fermions and two 
types of scalars, coupled in a symmetric way that employs two distinct scalar self-couplings 
and two Yukawa couplings, is analyzed. It turns out that it is possible to arrange for one of 
the Yukawa couplings to be asymptotically free. 

Dimensional regularization is employed to compute Feynman integrals for several 
reasons. First of all, massless tadpoles vanish when using this technique, eliminating a 
number of Feynman graphs. Secondly, divergences in an odd number of dimensions occur 
only beyond one-loop order. Thirdly, minimal subtraction is an efficient renormalization 
scheme which permits one to compute the renormalization group functions with relative ease; 
this technique is contingent upon using dimensional regularization. Fourthly, if one uses 
dimensional regularization in a three dimensional model that initially is free of dimensionful 
parameters, then one is not forced by consideration of radiatively induced divergences to 
introduce couplings involving massive parameters (such as masses for the scalar and spinor 
fields and quartic scalar couphngs). We now consider the various models. 



II. THE YUKAWA MODEL 

Initially, let us consider the model defined by the three-dimensional Euclidean action 

A 



S= d'x 



2 



(1) 



The conventions used are given in the appendix. 

In three dimensions, the couplings A and g are dimensionless. When dimensional reg- 
ularization im is used, only these couplings and the wave functions ip and (f) need to be 
renormalized. In computing the renormalization group functions, it is necessary to evalu- 
ate the divergent parts of the two-point functions and (■?/''?/'), the four-point function 
and the six-point functions (00^000). 

As was mentioned in the introduction, using dimensional regularization with this model 
means that tadpole and one-loop diagrams do not contribute to the renormalization group 
functions. Consequently, to order A and g^, we only consider the diagrams in fig.|l|, using the 
Feynman rules of fig.|^. In table | below, we list the number of distinct diagrams associated 
with each of the graphs in fig.|l|, the associated combinatoric factors, and the pole part of 
the diagram (in 3 — e dimensions). 

The divergences computed in table | do not contribute to vertices that occur in one-loop 
diagrams to the order which we are considering. 
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Analogous calculations are easily performed in the supersymmetric version of the model 
of eq.(|l|). The Minkowski space superfield action for this model is (using notation explained 
in the appendix) 



S 



$4 



(2) 



which in component form becomes 

S = J d^x [I [ADA + iij^'da^ijp + F'^ 



4! 



(3) 



where is a. Majorana two-component spinor in this theory; A and F are both real scalar 
fields. The equation of motion for the auxiliary field F may be used to re-express (§) in the 
form 



S = jd^'x \ [AuA + irdj'i^^ + ^AV' - 



(4) 



which is closely related to the action of (|1]). 

We determine the lowest order radiative corrections to the coupling A by computing the 
divergent parts of the integrals corresponding to the diagrams in figures |l|c, 0d, |l]e, p, and 
|I^ using the appropriate Feynman rules given in fig.^ The results are given in table |T|. 

We have also provided the results of computing the superfield diagrams of figs.|l|n and 
|l]o. We note that we could have examined the renormalization of the coupling A by com- 
puting the divergent part of the six-point diagrams rather than the four-point diagrams 
in fig.|l|; supersymmetry ensures that the results of this approach would yield an identical 
renormalization of A. 

A third model will now be considered. This will be done in Minkowski space and will 
incorporate two complex two-component spinor fields and '?/'2 which form the upper and 
lower components of a four component spinor field V', as well as a pair of real scalar fields 
A and B. The 2x2 Dirac matrices 7a in three dimensions, which are discussed in the 
appendix, are unitarily inequivalent to the matrices —7" (which of course satisfy the same 
Dirac algebra); the Dirac matrices we will use in this model are given by 



(5) 




= 0,1,2). 

We define two additional anticommuting Dirac matrices 

and Fs = 

The action which we will consider is initially expressed in the form 

S = J d^x[^l + Bp^B) - i^j^tp + f^^(A2 - B^ 

-tGi^T^MB - |(^' + B'') - ^ [a^B^ + A^B^ 





(6) 



(7) 
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where p = id and ip = ip~^To. There are two discrete symmetries present in this Lagrangian, 

B, B ^ -A, ip^Tr^ip, (8a) 

and 

A-^ B, B ^ A, ^ ^T^ij. (8b) 

These symmetries exclude couphngs of the form iptplA"^ + i?^), ijjipAB, i/jT^i/jIA'^ ± B^), 
'^V^ipAB, i^V^ilj^A^ ± ^2), A^B ± AB^, A^B^, A^B^ - A^B'^ and A^ - B^. They do not 
exclude '?/'r4r5^(^^ + B^), but this interaction is not perturbatively generated by radiative 
corrections induced by the action of (0). 

If in (^), we were to set G = g and a = A/5, then the action would become 



S 



d X 



-(f)*p^(j) - ^j/^ + gtfj [(f)^P+ + (0*)'P_) iJ - 



2^ 

90 



(9) 



where = (A+iB) / \/2 and P± = (l±r5)/2. The symmetry of (pa]) now becomes continuous: 



-ier5/2 



(10) 



as P±e-^^r5/2 ^ 

The interaction terms of (|^ are invariant under the transformations 



A= (A' + 5')/v2, B = [K-B' 



provided that we simultaneously redefine the couplings so that 

g' = G, G' = g, 
X — a 



a 



4 



, A' = (15(T + A)/4. 



(11) 

(12a) 
(12b) 



By exploiting the symmetries of ([l^), we can determine how G and a are renormalized to 
a given order in perturbation theory from results which give how g and A are renormalized. 
To do this directly we rewrite (0) in the form 



r I — ^ 

s= (fx -i<i>V$ - i^^^fi^i + \Y.^o. 

I a=l 



Q=l 



where we have reverted to using the two-component spinors ■j/'j and have defined 



A, B 
G, g ) 

150-+A 
4 ' 



A) 



(13) 



(13a) 
(13b) 
(13c) 



and 
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pa = ( 0-1, (73 ) . (13d) 

The transformations of equations (|rT|) and ( [T^ become 



i'^^{pi+p2)i^. (14a) 
$^^(pi + P2)$, (14b) 

^ /t2 Ci ^ C2 , (15) 



and 



under which the action of eg. ([T3|) is invariant. The Feynman rules associated with the action 
in eq.(|T3|) are given in fig-l^. 

In table |T| we provide the results of computing the relevant diagrams in fig.|l| in the 
context of the action of eq.(|l^). We have verified that these results can be obtained by 
using the action of eq.(0) directly, and also that the model of eq.(y) yields results that 
coincide with those following from eqs.(^ and (|l^) in the limit g = G, a = X/5. 

We are now in a position to use these results to determine the renormalization group 
functions associated with these models to lowest order. 



III. THE RENORMALIZATION GROUP FUNCTIONS 



The pole parts of the diagrams in fig.|l] as tabulated in tables |-|T| serve to fix the 
relationship between the renormalized and bare couplings, and hence the renormalization 
group functions, in the three models which we have considered In other words, we 

can determine the functions and in the equations 



9b = p' 



^R + l^ 

00 



9R + J2 



j/=i 



(16a) 
(16b) 



where, for concreteness we are using the model of eq.([l|). Similarly, the wave-function 
renormalizations are given by 



1 + E 



ct{^R,9R) 



(17a) 



and 

7 , , ^ ct{XR,9R) 

= 1 + 2^ — — • 

Altering the mass scale p so that 

p' = p{l + p) (p^O) 



(17b) 



(18) 
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leads to expressions for Xb and gs that contain contributions linear in e which can be 
eliminated by setting 

\r = Afl(l - 2ep) (19a) 

9r = 9Ril - ep) . (19b) 

If we now write and in terms of and gji and identify A'^ and g'^ (the renormalized 
couplings at scale /x') with terms independent of poles in e, we find that 

A'« = Xr - 2pa^ + 2pXR-^ + p~gR-^ , (20a) 

oXr ogR 

I ~ Q c \ daf _ dal , 
9r = 9r- pa-i + 2pAij-^ + PQrt^ , (20b) 

dXn ogn 

so that 

In determining to lowest order, one must keep in mind that a diagram involving a self- 
energy on an external propagator (so that the diagram is one particle reducible) has its 
divergence shared between the external wave-function and the coupling constant character- 
izing the vertex. This has the effect of reducing the contribution of the pole parts of the 
diagrams in figs.|I|f, |I^, and |I]1 to the appropriate coupling constant renormalizations by a 
factor of I in each of the models we are examining. 

Using table |, we find that for the Euclidean model of eq.(|I|), 

> (22a) 





2 

(22b) 

while for the supersymmetric model of eqs. (^) and (^ 

. 5A3 



487r' 



2 ' 



(23) 



as can be seen from table |T|. Finally, for the two component model of eqs. (|^) and (|T3|), 
table |TTT| gives the result 



af = 


9 


(3/- 


167r2 


af = 


G 


(3G2 


167r2 


= 


1 


f^A^ 
V6 


167r2 


< = 


1 


(3a' 


IGtt^ 



2G') , (24a) 
- 2g') , (24b) 



^a' + IQXg' + XG' + 15cxG' - 360g^ ] , (24c) 
i(TA + XG' + 15(tG' + 2Ag^ - 96G^ ) . (24d) 
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The results of eq.(0) are consistent with the symmetries of eq.(0). 

Together, eq.(pl]) and the functions Oi in eqs. (|2^^4D show that the rate of change of 
any couphng to lowest order is given by twice the corresponding ai, viz. for the Euclidean 
model 



dX 

dg 

for the supersymmetric model 



48 



45 



-9 

47^2 



(25a) 
(25b) 



d\ 5A^ 



and for the two-component model 

: 36-2 -2/ 
'5,o 5 



dg 


9 




87r2 


dG 


G 




87r2 


d\ 


1 




87r2 


da 


1 




87r2 



-A^ + -a^ + 16Xg^ + XG^ + 15aG^ 
6 2 



3a' 



\a\ + AG" 



360^^ 
15aG^ + 2Ag^ - 96G^ 



(26) 

(27a) 
(27b) 
(27c) 
(27d) 



We are now in a position to discuss the properties of these couplings. 
First of all, by (EsH), we see that 



9'o 



_9o_ In ( 

27r2 ^'^ I 



In/i/A 



or 9 = 



(28) 



where A is some renormalization invariant scale and g = go when /i = /io- Consequentially, 
the theory is asymptotically free, a property shared by non-Abelian gauge theory in four 
dimensions 0^ theory in six dimensions and (j)^ theory in four dimensions when the 
coupling is of the "wrong" sign [0. In the latter two cases, the models are unacceptable as 
they are energetically unstable. 

Since equations (|25|) are homogeneous in A, g^ and ln(/i/A)~^ we one can solve explicitly 
for the RG flow in the X-g plane and thus obtain a full solution to the scaling behaviour of 
both couplings at lowest order. In order to do this, we multiply (|25b|) by g and then divide 
by (^) to obtain the ODE 



dX 



5Ai _ ;^ 2 _ 45 4 
48 y 2 i' J 



(29) 



Then, by upon setting g'^ 
general solutions are 



Xz, it is straightforward to solve for -2(A). Unless g = 0, the 
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A = — -j= — for some k 7^ , or (30a) 



= Xz± , (30b) 



and 



where z± = (±y77/2 — l)/90. Notice that the generic solution ( ^Oa| ) appears to admit 
more than one phase; if k > then A varies smoothly with g and approaches zero in the 
ultraviolet limit. However if k < then A can branch outside of the perturbative regime 
when the denominator of ( |30a| ) is sufficiently small. Whether this occurs will depend on the 
particular choice of renormalization conditions. 

For the supersymmetric model, eq.(^) implies that 

l-T^ASln(^)=ln(A/.)- ^'^^ 

This model is not asymptotically free, in accordance with the general result of ref. P|. 

Finally, we consider the renormalization group equations for the two-component model 
in (13). By ( ^Tal) and (|7|) we see that 



This is a homogeneous equation; upon setting = zg'^ it can be solved easily to give either 



9' = G' 



47r2 



ln(A//i) 

(i.e. we have the Yukawa couplings of eq.(|^)) or 

g^'G^ <x g^-G^ 



(33a) 



(33b) 



Referring to figure ^ the graph of eg. (^3]), we see that the Yukawa sector has three distinct 
phases. \i G = g, then both couplings increase with scale according to (|33a|) with their 



equality preserved by the renormalization group fiow. \i g G then equation (|33b|) holds, 
and the smaller of the two couplings becomes asymptotically free while the larger evolves 
outside of the perturbative region with increasing momentum scale, in such a way that the 
combination + g'^ always increases with /i. To see this we set g^ = r sin 6 and G"^ = r cos 9 
in (|32|) to obtain 



dr 



3 — 5 sin 6 cos 6 
sin 6 — cos 6 



sin 9 + cos 9 
sin 9 cos 9 



(34) 



Another dominant feature of figure ^ occurs in the infrared region where the Yukawa 
couplings approach zero along the asymptote g = G. In fact, we have completed several 
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numerical solutions to the full four-parameter renormalization group flow of equations 
all solutions indicate that coupling configuration G = g, X = 5a (which leads to the chirally 
symmetric model of equations and (|10|)) is asymptotically realized in the infrared domain. 

In the special case when the Yukawa couplings g^ and are weak, then, in the ultraviolet 
limit, the six-point couplings tend to branch outside of the perturbative region parallel to 
one of the two UV-stable lines A = a or o" = 0. 

We can slightly alter the model of eq.(|^ so that there are Ng > 1 Fermions coupling 
with strength g to — B^, and Nq > 1 distinct Fermions coupling with strength G to AB. 
In this case the renormalization group equations of ( ^7a| ) and ( |27b| ) are replaced by 

{ANg + 5) - 2Ngg'^] (35a) 

{ANg + h)g^-2NGG^] . (35b) 



dG 


G 




247r2 


dg 


9 




247r2 



Again we find that there are special solutions in which the couplings flow in the following 
straight lines away from the origin as the energy scale increases: 

G = and g'^ = —— ^, , , , , , , (36a) 

{ANg + 5) ln(A//i) ^ ^ 

127r2 

g = and G^ = —— ^, _ , , . , , (36b) 

(4ArG + 5) ln(A//i) ' ^ ^ 



or 

G\fi) _ 6Ng + 5 
g^ifi) QNg + 5 



(36c) 



In addition there are generic solutions similar to the curves depicted in figure ^ in which 
either or G is asymptotically free depending on whether the fraction G'^/g'^ is larger or 
smaller than its critical value in 



IV. DISCUSSION 

We have considered a number of Yukawa couplings in three dimensions as well as the 
associated radiatively induced six-point scalar couplings. The lowest order contributions to 
the renormalization group functions for the couplings constants in these models have been 
computed, complementing the work of ref. [§, where a global SU{N) flavour symmetry is 
imposed on the interactions. 

One curious feature of these results is that the Yukawa coupling is asymptotically free 
in the model defined in Euclidean space by eq.(|I|), while the supersymmetric model in 
Minkowski space (with metric g^^ = (+ + ")) whose action is given in eqs.(|D and (^ 
is not. This apparent discrepancy is a consequence of the different properties of the Dirac 
matrices 7" in the two models. In Euclidean space the kinetic term ijjj/ip in eq.(|l|) is Hermi- 
tian provided 7" is identified with some unitary equivalent representation of the Pauli spin 
matrices a" (or —a"), in which case j/^ = ^^^q supersymmetric (Minkowski) model of 

eqs.(|]) and (|]), the kinetic term for the spinor field is Hermitian provided 7" is represented. 
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for example, by eg. (|A6|) , so that jf^ = —p^. This point of distinction is sufficient to alter the 
ultraviolet behaviour of the models, as the integrals associated with figs.|l]c-|l^ each contain 
spinor propagators. 

The two component model of eq.(^ has the peculiar feature that unless g = G, one 
of the couplings becomes small while the other grows as the renormalization group scale 
increases. This differs from the analogous situation in four dimensions where none of the 
couplings involving scalars and spinors is ever asymptotically free. 
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APPENDIX: CONVENTIONS 

In this appendix the notation and conventions that we have been using are outlined. 
In the Euclidean model, we work in Euclidean space with p = —id, and the Dirac gamma 
matrices are identified with the Pauli spin matrices so that 7^7* = 6""^ + ie°'^'^'j'^. This ensures 
that the kinetic term ip = ip^ . Feynman integrals are computed using 

(P^k ik'^Y 1 / ^yn/2+a-bVin/2 + a)Vih- a-n/2) , 
- 'm'^) -^—^ — (Al) 



(27r)" (F + m2)^ (47r)"/2 V J r{n/2)r{b) 

In the supersymmetric model, the conventions of are used. This means that an extra 

factor of (i) appears in eg. ([70]) , as the metric is given by g^^, = (+ H ). Indices are raised 

and lowered by the use of an antisymmetric tensor C^^, so that 

r = C''% , = ^^Cp^ (A2) 



where = -C°^ = M . Consequently, ii 5j = ^ ^ , then S^p = -Sj. A 



rank-two spinor Vap is identified with a vector {viz. equation (|A5|)) provided Vap = Vpa, or 
equivalently Va" = 0. Spinorial derivatives are defined by {9o,6"^} = {—ipa,6^} = ^J^ and 
the covariant derivative Da = da + iO^dap satisfies 



11 



□ 
1 



(A3) 



D'fiO) g{e) = d'x d'OfiO) D'giO) 



where = \D^D^ = -\D^D''. We also have 

These conventions are all consistent with replacing a vector VoP by 
where Va is the usual vector in 2 + 1 dimensions and 



i 
-i 



i 

1 



i 
-i 



(A4) 



(A5) 



(A6) 



with these conventions 



(7 



a 



(7 



a 



-ma 



13* 



and 



Wa\l')x^ + (7')/(7'^)a^ = -2^?"X^ . 



(A7) 



The wave equation for ip^ is now consistent with ip^ being real. The kinetic term 
ip°'{ida^)ipf3 = ip"'(j/)a^'ipi3 is Hermitiau. 

In the two-component model (|l5), some useful identities for the matrices pa in eq. (|13d|) 

are 



45 {pa®Pa ® a) [be 



and 



45 {pa ® Pa ® i; 



[bbbaaa] 



+ 12 {Pa)ba (Pa)aa ^ca + 12(p„)ca {Pa)aa ^ba 
+3 (Pa)a„ {pa)aa ^bc 



9 (Pa)6b (Pc.)b„ Saa + 9 (p„);,;, (pa)<,„ 5ba 

+ 18 (Pa)6a (P")6a '^^a + 9 iPa)ba (P«)aa '^^b • 



No summation is implied in these latter equations; square parentheses [■ ■ ■] indicate sym- 
metrization of all enclosed indices, while repeated Roman indices a and b indicate sym- 
metrization over specific indices (which are thus rendered indistinguishable). The coefficient 
45 is shown on the left hand sides to elucidate the 45 distinguishable permutations of the 
indices of p^ C?) p^ ® U given that both p„ and U are already symmetric. 
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FIGURES 






(k) (1) (m) (n) (o) 

FIG. 1. Feynman diagrams for the Euclidean, supersymmetric and two-component models 



-X 



-9 



-iDH{ei - 92) Ip^ 

-i\ 



FIG. 2. Feynman rules in the Euclidean and supersymmetric models 
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_a b -i(5"Vp2 

i ^ 1 

.a b, 

f ^^-.a-'' c 

e d 
a b 

i j 



FIG. 3. Feynman rules for the two-component model 




FIG. 4. Leading order flows for the Yukawa couplings in the two-component model of equa- 
tion (0) (The arrows represent the direction of increasing energy scale.) 
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TABLES 



Graph 


Number of diagrams 


Combinatoric factor 


Pole Part 


a 


1 


1 


^V/QOTT^e 


b 


1 


1 /O 

1/2 


g p /Wzir e 

3 /on 2 

— 5( /ozvr e 


c 


1 


1 


d 


1 


1/2 


3 //J /I 2 

— 5( /647r e 


e 


4 


1 




f 


2 


1/2 




g 


2 


1 


5A2/967r2e 


h 


10 


1/6 


i 


15 


1 


-ISA^VSSvr^e 


j 


90 


1 




k 


180 


1 


-455^/87r2e 
-A^VlBTT^e 


1 


6 


1 



TABLE I. Two-loop poles in the Euclidean model 



Pole Part 



-aVlBTT^e 




fA-V967r-e 
AV967 
-3iAV327r2e 
-aV487r2e 



TABLE II. Two-loop poles in the supersymmetric model 




Pole Part 



487r^e 



^(Atf-KlKi)(pi®Pl) + (l 



TABLE in. Two-loop poles in the two-component model. 



3^(Kf - Ki4){pi Pi) + (1^2) 
^(^f - ® Pi) + (1^2) 

^44^(76C? - 8C1C2 + 16C|)(pi pi 1) + (1 

Ti^(Ci«? + ^C2«i - i^CiKi)(pi PI u) + (1 

=^^{4 - nlnDipi PI U) + (1 ^ 2) 
+ KW2KP1 Pi(^ 1) + (1 ^ 2) 
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